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ABSTRACT 

Plasma turbulence is thought to be associated with various physical processes involved in solar flares, 
including magnetic reconncction, particle acceleration and transport. Using Ramaty High Energy 
Solar Spectroscopic Imager (RHESSI) observations and the X-ray visibility analysis, we determine 
the spatial and spectral distributions of energetic electrons for a flare (GOES M3.7 class, April 14, 
2002 23:55 UT), which was previously found to be consistent with a reconnection scenario. It is 
demonstrated that because of the high density plasma in the loop, electrons have to be continuously 
accelerated about the loop apex of length ~ 2 x 10 9 cm and width ~ 7 x 10 8 cm. Energy dependent 
transport of tens of keV electrons is observed to occur both along and across the guiding magnetic 
field of the loop. We show that the cross-field transport is consistent with the presence of magnetic 
turbulence in the loop, where electrons are accelerated, and estimate the magnitude of the field line 
diffusion coefficient for different phases of the flare. The energy density of magnetic fluctuations is 
calculated for given magnetic field correlation lengths and is larger than the energy density of the 
non-thermal electrons. The level of magnetic fluctuations peaks when the largest number of electrons 
is accelerated and is below detectability or absent at the decay phase. These hard X-ray observations 
provide the first observational evidence that magnetic turbulence governs the evolution of energetic 
electrons in a dense flaring loop and is suggestive of their turbulent acceleration. 
Subject headings: Sumflarcs - Sun: X-rays, gamma rays - Sumactivity - Sumparticlc emission 



1. INTRODUCTION 

In the standard flare scenario, the magnetic energy 
stored in the solar corona is released accelerating parti- 
cles. However, despite decades of effort, the exact mech- 
anisms remain elusive. It is still unclear how and where 
exactly the particles are energized. Hard X-ray (HXR) 
imag ing spacecrafts such as Yohkoh /HXT (|Kosugi et al.l 
19911 ) and RHESSI (|Lin et al.H200l are particularly use- 
ful for the spatial and s p ectral diagno s tics o f energetic 
parti c les dShibatal 11999): lAs chwandcn 2002 : iLin et al 



20031 iBrown fc Kontari 12005 
20101) . iMasuda et all (U99l . 
Liu et all (|2008f) : IJoshi et al 



ILin! 120061: iTemmer et a! 
199 % ISui et al.1 (120041) 
(|2009D have reported ob 



scrvations of a non-thermal HXR source above the top 
of an X-ray loop, suggestin g coronal energy relea se via 
magnetic reconnection (see iPriest fc Forbes! 120021 as a 
review) . 

Along with magnetic reconnection, another important 
element of the standard flare model is magnetohydro- 
dynamic (MHD) turbulence as it is believed to play 
a key role in a number of processes, from triggering 
fast reconnection to particle acceleration and transport. 
A considerable number of particle acceleration theories 
have been developed that require the presence of mag- 
netic fluctuations in the flaring loop. Assuming a high 
level of MHD waves, it has been shown that stochas- 
tic acceleration can effectively accelerate a large num- 
ber of both electrons and ions (|Miller fc Ramatvl 119871; 
Hamilton fc PetrosiaiJll992t iMelrosd 119941: iMiller et al l 



19971: IPetrosian fc Bykovl l2008t iBvkoy fc Fleishman! 
20091 iBian et alj|2010t iPetrosian fe Chenll2010D . These 



particle acceleration models include both resonant and 
non-resonant interaction between the particles and var- 



ious plasma waves. Waves and turbulence can be as- 
sociated eit her with current sh e ets d u ring magnetic re- 
conncction (IChiueh fc Zweibell U987t [S omov fc Kosugil 
119971: iLitvinenkoT 20061 ) or with reconnection outflows 
( Larosa et al.l Il994h . As the magnetic field in the solar 
corona is weak compared to the photosphere, direct mea- 
surement of even the mean field in the corona is prob- 
lematic. Therefore, the fluctuating components of the 
coronal magnetic fields and their relation to energetic 
particles remain unknown. In addition, the location and 
physics of turbulent acceleration within the general pic- 
ture of magnetic energy release is still awaiting an obser- 
vational basis. 

In this paper, using RHESSI X-ray observations and 
the visibility analysis technique, we show, for the first 
time, the presence of energy dependent cross-field trans- 
port of energetic electrons in a flaring loop. We estimate 
the relative magnitude and energy density of magnetic 
fluctuations which is required to explain these observa- 
tions. We show that these HXR observations provide a 
quantitative basis for the turbulence measurements and 
its governing role in the evolution of energetic electrons 
and are suggestive of turbulent electron acceleration in 
this loop. 

2. LIGHTCURVES, SPECTRA, AND SPATIAL 
DISTRIBUTION OF X-RAYS FROM A FLARING 
MAGNETIC LOOP 

We select a GOES M3.7 class flare, well-ob served by 
RHES SI , which was prev io usly an alyzed bv [Sui et al.l 
(l200l; iVeronig fc Browi] d200l : iBone et al.l (|2007l) ; 
IXu et all (|2008f) . ISui et"ai1 (|2004T > reported the existence 
of an X-ray source moving above the X-ray loop and con- 
cluded that the observations were consistent with the pic- 
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Fig. 1. — Time profile in different RHESSI energy ranges of the 
2002-04-14 event starting near 23:55 UT. The vertical lines indi- 
cates the three time periods for which the spectra, images and 
visibility forward fitting were performed. 



ture of the formation and development of a current sheet 
between the loop-top and coronal source. This April 14, 
2002 23 : 55 UT flare (Figures CD© shows a simple loop 
structure over the whole duration of the flare and, im- 
portantly, the magnetic loop is filled w ith dense plasma 
n ~ 10 u cm~ 3 (jVeronig fc Brownl 12001 ). Unlike the ma- 
jority of solar flares, which typically produce the soft 
X-ray (SXR) emission in the coronal part of the loop 
and bright HX R emission in the dense chromospheri c 
footpoints (e.g. lEmslie et alj|2003t iKontar et all 120081 ). 
the source of HXRs is coronal for this flare. Because of 
the high plasma density in the loop (Tabled]), energetic 
electrons of 10 — 20 kcV are collisionally stopped in the 
coro nal part of th e loop and produce thick-target emis- 
sion (|Brownlll971[ ). At the same time, for such a dense 
loop, 20 keV particles have collisional lifetime < 0.1s, so 
electrons need to be continuously accelerated in order to 
explain the observed HXR emission. 

We analyse this flare in three time intervals corre- 
sponding to the rise, peak and decay phases (see Fig 
[T|) . From the X-ray spectra (Fig [5J , we determine the 
electron acceleration rate dN/dt and the inverse power- 
law spectral index S of non-the rmal electron s, assuming a 
collisional thick-target model (|Brownlll971h . We also ob- 
tain from an i sothermal component in the spectral fit 
(jHolmanl [2005) : the electron temperature T e and the 
emission measure of thermal plasma, i.e. EM = n 2 V, 
where n is the plasma density and V is the volume oc- 
cupied by the thermal plasma. 

During the initial phase of the flare, acceleration of 
electrons up to ~ 50 keV is observed along with rapid 
heating of the plasma up to a temperature of 20 MK (Fig 
E}. As the flare progresses, the thermal plasma emission 
measure grows while the temperature slightly increases 
reaching 23 MK. The acceleration rate dN/dt is found to 
be around ~ 5x 10 35 electrons s _1 . In the decay phase of 
the flare, the thermal plasma starts to cool, the spectrum 
of energetic electrons steepens and the number of non- 
thermal electrons drops. The main flare parameters are 
summarized in Table [T] 

X-ray images reconstructed using the CLEAN algo- 
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Fig. 2. — The spatially integrated RHESSI X-ray spectra (black 
data points) for the three time intervals during the flare. These 
spectra were obtained us ing all front detectors excluding 2 and 
7 (sec Smith et al. 2002, for details) fitted (purple line) with 
an isothermal (red line) plus thick-target nonthermal (blue line) 
model. The pre-fiare background spectra are shown in grey. 



rithm (jHurford et al.l I2002T ) for the three time intervals 
are presented in Figure [3] A clear loop structure is ob- 
served up to ~ 20 keV. Since RHESSI measures the 
spatial distribution of X-rays only indirectly, via ro- 
tating modulating collimators (|Hurford et al.l 12002ft . X- 
ray images can suffer from their reconst ruction artifacts 
(|Prato et a l. 2009: Battaglia et al. 20TTI) and should not 
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TABLE 1 

HXR MAGNETIC LOOP AND FLARE PARAMETERS FOR THE EVENT STARTING APRIL 14, 2002 23:55UT. 



Time L(10keV)W(10keV) EM T e V (iTW 2 L/i)n = y/EM/V 6 dN/dt B^An/Bg 
(10 8 cm) (10 8 cm) (10 48 c m - 3 )(MK) (10 26 cm 3 ) (10 1Q cm- 3 ) (10 35 elec s" 1 ) (10 7 cm) 



23:57:48-00:04:00 
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21 
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9.3 
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2.3 


00:05:00-00:10:00 
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00:16:00-00:24:00 


23.5 


11.2 


6.9 


16 


8.9 


8.8 


9.8 


2.3 


$0.8 



be used to determine the spatial size of the source. 
Instead, a more direct and hence reliable approach is 
to 'stack' the rotationally-modulated lightcurves over a 
number of spacecraft spin periods to produce X-r ay vis- 
ibilities (jHurford et all 120021 ISchmahl et all 12001 . The 
latter are instrument-independent 2-dimensional Fourie r 
components of the HXR source (jSchmahl et al.ll2007h . 
For a simple source geometry, as displayed in this flare, 
the position and sizes can be determined unambiguously 
at various energies. An additional advantage of the X-ray 
visibility fitting approach is that the statistical errors can 
easily be propagated to find the uncertaint ies associated 
with the inferred parameters of the loop (|Kontar et al.1 
l2008tlPrato et al.ll2009t) . iKontar et~aTI (|2008ll2010D have 
successfully proved the robustness and sub-arcsecond 
precision of this method for HXR footpoints. X-ray 
visibilities are fitte d using a curved e l liptical gaussia n 
model of the flare (jHannah et al.l 120081 : Ku et al.ll2008h . 
This is done for a number of energy ranges between 9- 
20 keV and for the three time intervals (see Fig. [3]). 
The Full Width Half Maximum (FWHM) source sizes, 
specifically the length L(e), and width W(e) and their 
corresponding statistical uncertainties, are obtained and 
presented in Figure |U Using L and W at energy 10 keV, 
we find the volume V = ttW 2 L/A occupied by the ther- 
mal plasma. Assuming a volume filling factor of unity, 
the plasma density is calculated: n = y/EM/V (Ta- 
ble [IJ. It is f ound to be similar to t he previous den- 
sity estimates dVeronig fc Brownll20p: iBone et alj|2007t 
IXu et al.ll200l . 

As anticipated from the thick-ta rget scenario in a 
weakly changing plasma density (see I Brown et al.l 120021 
for details), the source length, L(e), in the direction par- 
allel to the guiding magnetic field, is a growing function 
of energy. The reason is that higher energy electrons can 
propagate further away from the region were they are ac- 
celerated, hence making the HXR source larger along the 
guiding field. In other words, in a dense coronal source, 
electrons travel a distance determined by their collisional 
losses. The stopping distance of an electron emitting X- 
rays of energy e while moving along the guiding field 
of the loop is ru ~ e 2 /(2Kn), where K = 27re 4 lnA, e 
is the electron char ge, In A ~ 20 is Coulomb logarithm 
(jBrown et al.l [2001 . For an energy independent length 
Lq of the accelerat ion region, the FWHM length of the 
source is given by (IXu et aUl200l : 



L(e) = L + a\\e 



(1) 



By fitting Equation ([TJ to the measured L(e), it is pos- 
sible to determine Lq and ct| [arcsec keV -2 ]. Figure @] 
shows that the length of the acceleration region Lq is 
around ~2x 10 9 cm at the peak of the flare. Moreover, 
since n — (2A"a||) _1 , au can be used as an independent 
measure to determine the plasma density. This approach 



gives plasma densities of 9.8 x 10 10 cm -3 , 2.9 x 10 11 cm~ 3 , 
1.7 x 10 11 cm~ 3 for the three time intervals analyzed in 
this paper. These densities are within a factor of two 
from the ones which were inferred above by using the 
thermal emission of the flare (see Table [lj . Larger den- 
sities at the peak and decay phases of the flare can be 
attributed to the presence of a colder plasma invisible in 
Soft X-rays (SXR) with RHESSI, as suggested in prev i- 
ous studies of thermal sources (e.g. Uakimiec et a"illl998D . 



3. ENERGY-DEPENDENT WIDTHS OF THE LOOP 

The CLEAN images (Figs. [3]) show that the loop width 
(size perpendicular to the guiding field) is like the length, 
growing with energy. This is seen quantitatively with the 
visibility forward fitted as a function of energy in Figure 
21 A growing width W(e) with energy is clearly present 
in the first and second time intervals (Figure[4j but could 
be absent in the third. This is particularly surprising as 
the guiding magnetic field is strong enough to ensure pure 
field-aligned transport of energetic electrons. Indeed, for 
a mean magnetic field of ~ 150 Gauss as measured in this 
loop using radio data (jBone et al.ll2007l ). the Larmor ra- 
dius of tens of keV electrons is ~ 2 cm and, therefore, 
the trajectories of their guiding centers coincide with the 
magnetic field lines. Collisions are unimportant in pro- 
ducing perpendicular transport for the parameters of the 
plasma, only resulting in a radial excursion (classical col- 
lisional diffusion) of the order of the electron gyroradius 
in a collisional time. 

However, it is well-known that particle transport across 
a mean magnetic field can be dramatically enhanced, well 
above the collisional value, in the presence of turbulence. 
The reason is the braiding of the field lines themselves 
due to perpendicular magnetic fluctuations. In the pres- 
ence of magnetic fluctuations perpendicular to the 
mean field (Bo = Bqz), perpendicular transport of field 
lines occurs which is governed b y the following equation: 
dv ( z)/dz = B (r i (z) , z) /Bp (jRechester fc Rosenbluthl 
1978). Field lines braiding is believed to be a dominant 
proces s governing t he transport of particles i n the solar 
wind (|Jokipiil Il966l ). in laboratory plas mas dBickerton 
1997 ). and in thermal coronal loops dGallowav et al ' 
2006). It is safe to neglect the contribution of per 



pcndicular electric fluctuations and subsequent E x B 
drift to cross-field transport. For MHD fluctuations, the 
electric field is E± ~ vaB±/c such that the magnetic 
drift vg = vu(Bx/Bq) dominates over the electric drift 
ve = cF,± x B /i?o ( c is the speed of light and va is 
the Alfvcn velocity). Since vb/ve ~ v \\/va, the contri- 
bution of the E x B drift to perpendicular transport is 
negligible for non-thermal electrons with energies above 
> 10 keV (uii ~ 6 x 10 9 cm s _1 ), i.e. Vu 3> va, with va — 
10 8 cm s" 1 in the magnetic loop. The magnetic field 
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Fig. 3. — CLEAN imago ovcrlayed by 30, 60 and 90% levels from the forward fitting to X-ray visibilities. Each row indicates a different 
time range. Each column a different energy range. The parameters of the loop are summarized in Table 1 and Figure [4] 



line diffusion coefficient is given by the Taylor formula: 
D M = {l/Bl)^dC < B ± (r ± (C),C)Bi(r ± (0),0) >. 
Here, < Bj_(rj_, z)Bj_(0, 0) > is the two-points corre- 
lation function and <> denotes an ensemble average. In 
the quasi-linear approximation, the perpendicular trans- 
port of field lines is characterized by a diffusion coefficient 
Dm — (Bj_/Bq)\u, where Am is the parallel correlation 
length of the perturbations B j_ (jJokipii fc Parker! 1 19691 : 
iRechester fc Rosenbluthlll978l ) . An electron moves a dis- 
tance r\\ ~ e 2 /2Kn along the mean field in a collisional 
life-time. Therefore, this electron also makes a radial 
excursion given by r± = y/2D~M~r~\\ , as a consequence of 
the radial transport of field lines. This means that the 



FWHM width of the source W(e) should also grow with 
energy. It is the sum of the acceleration region width Wq 
and the diffusion part due to perpendicular transport: 

W(e) = Wo + cue, (2) 

where a± = ^JTD~mo\\ is measured in arcsec keV -1 . The 
unique feature of this expression is that it depends on 
the level of magnetic fluctuations via the magnetic field 
line diffusion coefficient Dm- Therefore, by measuring 
the energy dependent width of the flaring magnetic loop 
(Figure 4), we find automatically the diffusion coefficient 
for the magnetic field lines, i.e. Dm = \\\B\_/Bq = 
a\/(2a\\) (Table 1). As anticipated, W(e) can be well 
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fitted with a linear fit (Figure HJ to find a± and the 
magnitude of magnetic fluctuations. The magnetic dif- 
fusion coefficient \hB±/Bq is changing with time and 
is strongest when the acceleration rate dN/dt of non- 
thermal electrons is at the maximum and the spectral in- 
dex of the accelerated electrons is smaller (hardest HXR 
spectrum). The magnitude of fluctuations at the decay 
phase is weaker, so the data cannot provide a reliable 
estimate, but only an upper limit. 

The energy density of magnetic fluctuations in the loop 
is B\ /8tt ~ {B%/8it)D m /\\\. Taking the length of the 
acceleration region 2 x 10 9 cm as the upper limit for 
correlation length Am, the energy density of magnetic 

fluctuations is ~ 10 erg cm -3 and the relative level of 
fluctuations Bj_/B ~ 0.1. For Ay = 2.3 x 10 7 cm, the 
magnetic perturbations follow to be strong B±/B = 1, 



and the energy density of turbulence ~ 900 erg cm -3 . 
This value is close to the energy density of the flaring 
plasma 3nksT e ~ 830 erg cm -3 and much larger than 
the energy density of the non-thermal electrons 



U, 



Hi 



F(E) 



-AE, 



(3) 



where S = ttWq/4 is the cross-sectional area of the loop, 
F(E > E c ) oc E~ 5 is the thick-target acceleration rate 
differential in energy [electrons s _1 keV and E and 
m e are the electron energy and mass respectively. Using 
the values from thick-target fit and loop width measure- 
ments, the energy density of non-thermal electrons, U n th 
is about 11, 6, and 1 erg cm -3 for the three time inter- 
vals respectively. Therefore, we can conclude that the 
energy density of the magnetic turbulence is sufficient to 
energize at least non-thermal electrons and possibly the 
whole flare. 

4. DISCUSSION AND CONCLUSIONS 

X-ray visibility based analysis has allowed us to infer 
the energy dependent source size for a well observed flar- 
ing loop. Both the length and the width of the loop are 
found to be larger at higher energies. This is in clear 
contrast with HXR footpoints, which show the opposite 
trend; higher en ergy sources appear smaller than lower 
energy sources dKontar et al.l [20081: iPrato et all 120091: 
iKontar et al.ll2010l) . An energy dependent growth of the 
loop length is expected from the colli sional precipitation 
of energetic ele ctrons along field lines (jBrown et alj|200"2l : 
IXu et aiT 200cf) . The high density of the loop implies that 
the electrons are accelerated inside the loop region with 
the characteristic size 2 x 10 9 cm. At the same time, 
the size of the loop perpendicular to the strong magnetic 
field is an increasing function of energy. The increase 
of the loop width implies strong cross-field transport of 
non-thermal electrons and is found to be consistent with 
the diffusion of magnetic field lines. This allows us to 
determine the level and locations of magnetic fluctua- 
tions in the flaring loop. For comparison, we note that 
the magnetic diffusion coefficient Dm = X\\Bj_/Bq mea- 
sured in situ in the solar wind, where B±/Bo ~ 1 and 
Am ~ 10 11 cm (Matthaeus et a.1.1 11986T ) is about four or- 
ders of magnitude larger than our values presented in 
Table CD 

We note that although the length and width of the 
flare are likely to be governed by parallel transport and 
diffusion of the magnetic field lines, equations (|H2j) for 
the width and the length are approximations. The si- 
multaneous treatment of acceleration, energy losses and 
cross-field transport is needed to assess the particle dy- 
namics more accurately. This might lead to super or sub- 
diffusive cross-field transport of electrons and cannot be 
ruled out by the current data, so the radial excursions 
can scale as r± oc e' 3 , with /3 close to one. Nevertheless, 
the width-energy measurements provide a powerful and 
unique tool to measure magnetic fluctuations in a flar- 
ing plasma. Importantly, the energy density of magnetic 
fluctuations appears to be energetically significant and is 
found to be larger than the energy density of non-thermal 
particles. 

The temporal and spatial analysis by iSui et all ([2004D 
of the weaker source above the dense loop analyzed here 
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suggests a picture consistent with formation of a recon- 
necting current sheet between the loop top and the coro- 
nal source. Assuming this reconnection scenario, our ob- 
servations suggest the presence of MHD turbulence in 
the loop under the reconnection site. During the event, 
the length of acceleration site increases only by ~ 10%, 
while the width of the bright HXR loop Wo grows with 
time more than a factor of ~ 2, which could be attributed 
to continuous piling up of newly reconnected lines. The 
plasma mode responsible for the perturbations remain 
unknown, but these observations indicate the presence 
of noticeable magnetic fluctuations during the solar flare 
inside the loop, where the acceleration of electrons is hap- 



pening. It is possible to speculate that the magnetic 
field fluctuations are produced by magnetic reconnection 
above the loop and accelerate electrons inside the flaring 
loop. 
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